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We construct the effective transition operators relevant for neutrinoless double electron capture 
leading to final nuclear states different than 0"*". From the structure of these operators we see that, 
if such a process is observed experimentally, it will be very helpful in singling out the very important 
light neutrino mass contribution from the other lepton violating mechanisms. 
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Csl . INTRODUCTION 

^ ■ The nuclear double beta decay is process considered very long time ago[l|. It can occur whenever 

■^r [ the ordinary (single) beta decay is forbidden due to energy conservation or greatly suppressed due 

to angular momentum mismatch. It could proceed with the emission of two neutrinos or in a 

neutrinoless mode. The 2^'/3/3— decay is just an example of a second order weak interaction. The 

i~~i' exotic neutrinoless double beta decay (Oi^/3/3— decay) is the most interesting since it violates lepton 

'^ [ number by two units. It was first considered by Furry [3| more than half a century ago as soon it was 

II ^ realized that the neutrino might be a Majorana particle. To-day, seventy years later, Oi^/3/3-decay: 

5^. {A,Z)^ {A,Z + 2) + e- +e-, (Oi^ /3/3 decay), (1) 

A = Q, continues to be one of the most interesting processes, since: 

^ . •It will establish whether there exist lepton violating processes in the Universe. 

\o ' . 

Qv ' • It is perhaps the only process, which can decide whether the neutrinos are Majorana particles, 

^S| . i.e. the particle coincides with its charge conjugate, or Dirac particles, if it does not. 
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• It will aid in settling the outstanding question of the absolute scale of the neutrino masses, 
especially if this scale turns out to be very small. Neutrino oscillation experiments can only 
determine the mass squared differences. 

For recent reviews see [3j-iQ] . 

If this process is allowed other related processes in which the charge of the nucleus is decreased 
by two units may also occur, if they happen to be allowed by energy and angular momentum 



;h ' conservation laws, e.g. 

{A,Z)^ (A, Z - 2) + e+ + e+ (Oi/ positron emission), (2) 

A = Q -AnieC^, 

{A,Z) + e^ — > {A,Z — 2) + e^ (Oi^ electron positron conversion), (3) 

A = Q -2meC^, 

{A, Z) + e^ + e^ — > {A, Z — 2) + .xrays (Oz^ double electron capture ), (4) 

A = Q, where A is the available energy. 
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The last process is always possible, whenever ([3]) is. It has been first considered long time ago 
[71, y as a two step process: In the first step the two neutral atoms, (A,Z) and the excited (A,Z-2), 
get admixed via the lepton number violating interaction . In the second step the (A,Z-2) atom 
de-excites emitting two hard x-rays and the nucleus, if it is found in an excited state, de-excites 
emitting 7-rays. The life time expected was very long, since the above mixing amplitude was tiny 
compared to the energy difference of the two atoms involved. It has recently, however, been gaining 
in importance |9l4l2[ after ion Penning traps [l3| made it possible to accurately determine the Q 



values, which gave rise to the the presence of resonances, which in turn could lead to an increase of 
the width my many orders of magnitude [9|. Decays to excited states are in some cases possible and 
provide additional experimental information, e.g 7 rays following their de-excitation. 

Furthermore neutrinoless double electron capture is of experimental interest, since one expects the 
reaction to be essentially free of the two neutrino background, which plagues ordinary neutrinoless 
double beta decay. It has, also recently been realized [10| that the resonance condition can be 
fulfilled for transitions to excited states. Thus, e.g, in the case of ^^^Dy— >^^^Gd, such excited states 
canbe^ 0-,l~,0+,l+ and 2+. 

Before proceeding further in examining neutrinoless double electron capture, we note the following: 

• Transitions to 0~, 1~ states have not been considered previously in standard double beta decay, 
since they are energetically forbidden. Such transitions can naturally occur, however, in the 
presence of right handed currents, arising by interference between the left handed and right 
handed currents, i.e. of the type jl x Jr, with an amplitude is essentially independent of 
the neutrino mass. We will show that this process can occur in double electron capture, even 
though the electrons are bound. 

• Transitions to 0^ states, both the ground states and excited states, can occur both in the mass 
mechanism and the jl x jr. The latter contribution in ordinary double beta decay occurs since 
the produced electrons have sufficiently high momenta. In double electron capture instead of 
the electron momentum one can exploit the nuclear recoil term. 

• transitions to 1+ and 2+ can occur only, if the electron momentum is sufficiently high, as in 
neutrinoless double beta decay, or by retaining the nuclear recoil term in the hadronic current. 
In double electron capture this can naturally happen, if one of the electrons is captured from 
a p-orbit. 

In the present paper we will derive the effective transition operators relevant for neutrinoless double 
electron capture, emphasizing the mode which is independent of the neutrino mass. The neutrino 
mass contribution has already previously examined, see e.g. the recent work and references therein 
i- 



NUCLEAR TRANSITION OPERATORS LEADING TO NEGATIVE PARITY STATES 

To leading order such transitions cannot occur in neutrino mass mechanism. In the j^ x j^ case, 
however, the chiralities involved are such that in the intermediate neutrino propagator one picks the 
momentum, not the mass. The neutrino momentum leads to the following: 

• The time component of the momentum leads to an average energy difference -En.i ~-£'„,2, where 

En,i =< En > Mi + Ee{i), with Ee{i) the electron energy, i = 1,2. This can lead to positive 

parity transitions in double beta decay, but it is not the dominant contribution. 

• the space part of the momentum leads to an effective operator which proportional to f^- 
with r = ri — r2, r = |ri — r2|. This naturally leads to negative parity transitions. 



^ After completion of this work, two papers have been published Il4l uM , which show that the resonance condition 
can be met for such excited states 



Positive parity transitions can, of course, occur in neutrinoless double beta decay via this mechanism 
by retaining in the dipole term in the expansion of the essentially distorted plane wave function of 
the two electrons, which leads to the term i (pi.ri + P2-r2) '^12^7- This possibility does not exist in 
neutrinoless double electron capture. We will, however, examine alternate possibilities in the next 
section. 

The effective lepton violating transition operator, within the closure approximation over the in- 
termediate states, which holds for neutrinoless double beta decay and is expected to be even better 
in neutrinoless double electron capture, can be cast in the form: 

^ ^ ^ ^ ^ n,m 

with (/'ni.o(^n)0n2.o(^m) (^n = ^1 j ^m ~ ^2) thc bound s-clcctron wave functions and 

ujl,,k = u(pi)7M7fc7i-(l - I5)u''{p2), u^^fc = u(pi)7^7fe7^(l + 75)u''(p2) (6) 

^lH = T-(n){gv, -5ACT„), J^{m) = T_(m)(.gy, +5^^™) (7) 

/(a;„,m) = ^--J((Jo), X ^ x{n,m) ^\rn-rm\/Ro (8) 

ax X 

2 l"^ sin ?; 
JiSo) = - / ^-rfy, <5o = (^ ^„ ^ -M,)\r^ - r„| (9) 



in double electron capture: 



(5o ~ 0, J((5o) « 1 and /(a;) « -^ 



For non relativistic electrons one finds that the effective transition operator can be cast in the 
form: 

f . An operator associated with the time structure of the leptonic current: 

4 = w(pi)7o75u''(P2), (10) 

which is 

^0 = 7rm—fy}.^9\ y^ (t)mfi{rn)4in^fi{rra)T-{n)T-{m)fnm-'jJo{n, rn)f{xn,m), uJoin,m) = J(cr„XCTm) 

^ ^' n.m 

(11) 

which can cause 0+ — > transitions. 
2. An operator associated with the space structure of the leptonic current: 

/= u{pi)ju'{p2) (12) 

that is 



G? 



I 



E'?^"i.o(^«)0n2,o(»'m)2T_(n)T_(m) [f„„ X uji(n,m)] /(x„^„i), uJi{n,m) = (cr„ - cr„) 

(13) 



Since the Os electron wave functions vary slowly inside the nucleus one can replace them by there 
average value and write: 

fl, = A.« ^ (Os)2a| ^ a, A„,.= ^^^^2^3 (14) 

where the A„ff has dimensions of energy and the average is over the radial part alone (the integral 
over the angles is unity). The dimensionless effective transition operators now become: 

^0 = 45^^r_(n)r_(TO)f„m.a;i(n,m)/(a;„,m), uJo{n,m) = i{an x am) (15) 

n,rn 



Ql =2gvgA^T^{n)T^{Tn)fnm'X(^l{n,Tn)f{Xn,m), LL>i{n,m) = an- (Tm (16) 

n,rn 

One can show that averaging over the initial polarizations and summing over the final ones one gets 

1 

|M^|2 = 2^1 ^ /||r2f||0+ ^ P, fc = i = J/ = rank of the operator (17) 

4=0 

where the factor of 2 comes from the leptonic current, k = Jf is the rank of the operator and the 
double bar indicates its standard reduced ME, which can cause 0+ — > 1^ transitions. The mixing 
V"^ [ij] between the two atoms takes the form: 

V^ = Al, (^ {Osfa% yf \ME\^ (18) 

TRANSITIONS TO POSITIVE PARITY STATES IN ELECTRON CAPTURE 

We will distinguish two cases: 

The captured electrons are in s states 

As we have mentioned the 0+ -> 0+ transitions associated with double electron capture have 
already been studied during the last few years. The contribution arising from the neutrino mass 
independent mechanism has not been studied since the electron wave function cannot provide the 
extra parity change. This , however, can be accomplished via the recoil term in the nucleon current^. 
In other words one of the hadronic currents is of the form discussed above, but the other is of the 
form [16|: 

J^in) = T^{n){gv'Dn, -gACn), JR{m) = T^{m){gv'D.ni, +gACrn) (19) 

with 

C„ == o-„.(g„- 2F„), D = - [{qn - 2Pn) ^ ian X Qn] (20) 

Zttin IniN 

Where rriN and P„ are the nucleon mass and momentum of nucleon n, while g„ is the momentum 
transfer to that nucleon. 

The procedure is a bit tedious but one can proceed as above to arrive at 



^ It is interesting to note that in the case of neutrinoless double beta decay for the mass independent contribution 
one needs a correction either from the dipole term of the electron wave function or the nucleon recoil. Both are 
first order corrections. The nucleon recoil term has not been extensively studied even though it has been suggested 
by Kotani al long time ago [16| | 



1. One operator associated with the leptonic current: 

f=7l(pi)775U^(p2) (21) 

which is associated with an operator appropriate for 0+ — > 1+ transitions, namely: 

^r"' = 7^-j;:^2^0„i,o(''n)0n2,o(?'m)T-(n)T_(m)(a;„™ +a;;„)/(a;(n,m)) (22) 

n^7n 

with 

W„m = gvgArnmiCn — Cm) (23) 

W'mn = ^^^(Dn - Dm) X fn,m (24) 

2. One operator associated with the leptonic current: 

£5 = u{pi)^5u''{p2) (25) 

which is suitable for for 0+ — > 0+ transitions, namely: 

^r°" = ^:5;:^2 2j'/'«i,o(r„)0n2,o(?-m)T-(n)T_(m)(a;,„„ +a;;„) (26) 

^7un = gv9Arnm- [«(-0-„ X D,„ + (T„ X £)„)] (27) 

and 

^nm = 9A^nm- [—C^nCm + C^mCn] (28) 

Under the factorization approximation one can write: 

nr°" = A„„ -< {Osfal >- nr°'\ i = 0,l (29) 

j^„co.i = 2^T_(n)r_(m)(c^„„ +cj;„)/(a;(n,m)), (30) 

^rccoa = 2^r_(n)r_(TO)(w„„ + w;,„)/(a;(n,TO)) (31) 

and then use Eqs ^ and (fTSi) . 

There remains the problem of transforming the operators C„ and D„ into coordinate space. The 
operator q„ is local and leads to the operator —i^^-^^ where x = t/Rq, r being the relative distance 
of the two nucleons. The operator P„ is non local. Thus one finds: 

Pn -^ -*;^ (v„ - V„) (32) 

where the operator v„ operates with respect to x„ on the nuclear wave functions on the bra and 
ket as indicated by the arrows. As a simple example we will rewrite the operator the operator 
fnm{Cn — Cm) in Coordinate space. We find: 

^nva\S^n ^m) 77TT T^^nni i ynm-^n ^7im-'^ni)~j 

lAinKo I ax 



{Vnfi^) - /(a;) Vn) ■<^n - {ymfi^) " f i^)^ m) ■'^rn } 



(33) 



This can be rewritten by showing the tensor character of the various components as: 
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2^2 / 



Vs V 



^yyjji ^ jri ) r— I [' 'tun "• • rirn] ^ \^n ^m) I 

(v„/(a;) - /(a;) V„) o-n - J'nm- (Vm/Ca^) - /(a;) Vm) f^^ 
Ti ® ( V„/(a;) - /(a;) V„) 
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rnm^" Vm/(a;) - I{x)SJ 



In the case of the local scalar spatial contribution it is understood that: 



dl 

dx 



Attx"^ 



(34) 



(35) 



The overall tensor rank is J = 1, the spin rank is one, but it contains both scalar and tensor 
components in the orbital sector. It can thus lead to 0+ — )• l"*" transitions. The decomposition of the 
other operators, if needed, can be done in a similar fashion. It is quite tedious but straightforward. 
The scale of the nucleon recoil contribution is set by (l/2TOjvi?o) = 0.07A~^^^. In the case of 
neurinoless double beta decay the square of this scale should be compared with the ratio of the 
phase space of this jl (X" Jr going through the electron momenta to that of the mass term. 



One of the electrons is in a p orbit. 

If the resonance condition is met for such an electron, one may consider this possibility, even though 
the probability for finding a p-electron inside the nucleus is tiny. The factorization approximation 
should not be used in this case and the electron wave function must be part of the nuclear transition 
operator. The average over the nucleus of the square of the electron wave function is, however, still 
a good measure of the scale of the expected nuclear matrix elements. Thus combining the above 
operator for the decays to negative parity states we find that one can have transitions to 0+, 1+ and 
2+. The spin structure of the operator is still: 



CJo = lO-n 



' dm or UJi = o-„ 



(36) 



i.e. antisymmetric of rank one. There appears to exist a wealth of orbital contributions. 

Indeed let us assume that the p electron has total angular momentum j and the initial angular 
momentum of the atom is J. The structure of effective transition operator can be cast in the form: 



ioo-fio ^ C/(l, 1/2, i, 0, 1/2, 1/2, 1, 0, J)Sj^i 
for the time component and 



uJo.rV 



4 



(37) 



(wi. X f).e^ C/(l,l/2,j,0, 1/2,1/2,1,0, J)L/(1, 1,0, 1,1, J, J, 0, J) 
{lui X f)(E)V e.s 



(38) 



for the space component {£ is the leptonic current not to be confused with the orbital angular 
momentum operator). U{...) are the usual unitary nine-j symbols needed for the re-coupling of the 
angular momenta involved and V essentially is the two electron orbital wave function: 



V 



1 V3 1 



47r a/Stt V2 ^ 



(a" 



X2 + A Xl 



(39) 



with 

p p 

A ^ —{(l>ni,o{rn),(l>n2,l{rm)), A = — (0„j j (r„), 0„2 ,o(?'m)) (40) 

with (j)ni,ei, i = 1, 2 the radial electron wave function. We have chosen to express the length in units 
of the nuclear radius Rq, i.e. x — r/ Rq, with compensating factors in the definition of A and A , 
a procedure quite standard in the evaluation of the nuclear ME. 

We found it convenient to go into the relative x = Xn — Xm and center of mass 2X = x„ + Xm 
coordinates and thus we get: 

V^ — (A(n,m)x + B(n,m)X), A(n,m) = Vs^- -, B(n,m) ^ Vs^iA' + A) 

An ^ ^ ' y ^ ' " \ 1 / ^ 2 \ ^ / ^\ i 

(41) 

The operator Y must be symmetric under the interchange of particles n and m to yield a total 

operator which is overall symmetric under the exchange of the particles n and m. So only the 

B-terms survive. 

The next step is to perform an additional re-coupling bringing together the two spacial parts. One 

then finds: 

1. For the time component: 

,.^4^ :^(^^Co(j, L)[^o®r^(r,i?)]^<5j,i4 (42) 



4. 



with 



and 

T^=r.y., T^=fxX, r2 = [f(^X]2 (43) 

2. For the space component: 

^1 X f^o -^ ^^^^ Yl ^2 U, J, L)[uJi<E> T^ (f , R)] ^Is (44) 

L 

where the coefficients C2(i, J, i) are given in table H] 

By combining Eqs (fTTj) and (J42|) as well as Eqs (fT3|) and (|44l) we obtain the complete operators 
relevant for the decay of 0+ to positive parity states with J = 0, 1 and 2. 

To simplify matters one may seek an approximation familiar from the Os case, i.e. to take the 
radial part of the electron wave function outside of the integrals and replace it with its average value 
over the nucleus. Thus: 

~(2 



a=A„„y^4B^^„ A.„. = 5^ J 3 , BoB(n,m) (45) 

' 2v2 "'^o'^B 

where a^ is Bohr radius and, again, the averaging is over the radial part alone. Thus the effective 
transition operators now become: 

^0 =4g^^T_(n)T_(m)/(x„,„)^Co(j,L)N®r^(r,i?)]''<5j,i (46) 

n,m L 

f^i =2gy5^^r_(n)T_(m)/(a;„,™)^C2(i, J,£)[wi®T^(f,i?)]^ (47) 



In fact it will be adequate to consider 



-,2 . , , /r, ^2 , , r, I Rq 



r 

where Rn.e are the radial parts of the bound electron wave functions. Note that the p wave function 
involved in the averaging has been stripped of an extra power of r, since the r has been incorporated 
into the nuclear matrix element. 

The expression for the nuclear ME is the same with that of Eq. (|17p , where k = J = Jf is the 
rank of the operator and the coupling of the two atoms is given by: 

^2 ^ A^,.^ ^ a|sV \ME\'^ (48) 

We will check the factorization approximation using realistic atomic 3s and 4^3/2 wave functions 
provided to us by Shabaev [l7|, which were obtained by standard methods (see, e.g., Grant [l8| and 
Bratsev et alfi^). 

• Using the density profile shown in Fig. [TJ corresponding to a phenomenological normalized to 
unity Woods-Saxon density, we find: 

^ a|j(i?3,) V= 2.4 X 10^ ^ 3a| U— i?4p3/ J )-= 3.2 x lO^^ — y 

2 



-< 3 (^Rip,^,j >- I < (Rssf ^= 1.3 X 10- 



• Using the density profile corresponding to a simplified S.M. wave function shown in Fig. [2l we 
find: 

^ a|j (i?3,) V= 1.4 X 10^ ^3al(—Rip^A ^= 1.8 x 10"^ — > 

^ 3 {{^R^pA ^ / -< (i?3,)2 y^ 1.4 X 10^^ 

Had we not used the factor 3{Ro/r)'^ the p-integrals would have been quite a bit smaller. The 
agreement in the ratio of 4p to 3s probabilities between the two nuclear densities is quite good. Fur- 
thermore, using a simple Os HO wave function, we find that the factorization approximation works 
to within 10%. More detailed computations using realistic and more accurate shell model wave 
functions are currently underway. The results, with and without employng the factorization approx- 
imation, will be published elsewhere [20|. The nuclear matrix elements leading to the ^^^Dy— >-^^^Gd 
leading to sharp resonaces 15] 0~, 1~, 0+, 1+ and 2+ will also be calculated ^], using realistic wave 
functions both for the intial and the final nuclear states. 



DISCUSSION 

As we have mentioned the observation of neutrinoless double decay will be a great event in 
particle and nuclear physics. It will establish the existence of flavor violating interactions and it will 
demonstrate that, independently of the dominant mechanism causing it, the neutrinos are Majorana 
particles. More importantly, it is expected to complement the great successes on neutrino oscillations 
by determining the elusive up to now absolute scale of neutrino mass. For this last great contribution 
to become unambiguous, it will be necessary to determine what fraction of this decay is due to the 
light neutrino mass compared to all the other competing mechanisms ^]. That is a very difficult 
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Figure 1: The (normalized) Woods-Saxon density profile employed in obtaining the average 

electron wave functions. 




Figure 2: The nuclear S.M. density profile employed in in obtaining the average electron wave 
functions and checking the factorization approximation. 



task indeed. It may become feasible, if there exist data in as many targets as possible. In fact it has 
recently been shown that, within the mass terms alone, even if both left and right handed currents 
are involved, this is, in principle, possible [2ll . [22 1. Provided, of course, that experimental data 
involving a sufficient number of targets become available. It will be much more difficult to achieve 
this, if the mass independent contribution involving the jl x jji contribution becomes important. 
In this direction the observation of double electron capture will be very useful. More specifically we 
have seen that: 



Table I: The coefficients C2(j, J, L) discussed in the text 



C2iJ,J,L) 

^T 

2 

_J_ 
2^6 



i 

4 
_2 



1. The neutrino mass term can proceed in neutrinoless double electron capture as in ordinary 
double beta decay. For the latter to be practical, of course, the resonance condition must be 
met. 

2. Neutrinoless transitions to 0+ ^ J^ 7^ 0+ can occur only if there exist right handed interac- 
tions through the right left interference in the leptonic sector (jl §5 Jr)- Then: 

• The unusual transitions to the negative parity 0^ and 1^ sates can naturally occur in 
double electron capture, since they are energetically allowed, if the resonance condition is 
met. So such an observation will unambiguously determine the neutrino mass independent 
lepton violating parameter, which is hard to get in neutrinoless double beta decay. 

• Transitions to positive parity states 7^ 0+ in ordinary double beta decay can occur either 
through the dipole component of the electron wave functions, the most favored scenario, 
or if one of the hadronic currents involved contains the parity changing nucleon recoil 
terms. 

• Transitions to positive parity states with angular momentum less than 3 in neutrinoless 
double electron capture can proceed if a) one of the hadronic currents involves the recoil 
term, while both electrons are in s-statcs or b) more naturally via the negative parity 
operator involving ordinary nucleon currents, if one s-electron and one p-electron are 
involved in the capture. The probability of finding a p-electron in the nucleus is, however, 
quite small, but this disadvantage maybe overcome by the collaborative effect of the 
expected large nuclear matrix elements and the resonance condition for capture to 2+. 

• The transition operators involved have a reasonable structure, but nuclear matrix ele- 
ments of such operators have not yet been calculated for any system. We have no reason 
to expect them to be suppressed. The nuclear structure calculations involving such nu- 
clear systems may be a true challenge to nuclear theorists. 

It is clear from the above that the neutrinoless double electron capture can distinguish between the 
mass term contribution and the neutrino mass independent jl x jji interference term. It can do 
it more efficiently than ordinary beta decay due to the possibility of the negative parity nuclear 
transitions. It cannot, of course, distinguish between the various mass terms (light neutrino, heavy 
neutrino, SUSY contribution etc), but neither can ordinary double beta decay. This can only be 
done by exploiting the results of at least four different experiments (leading to 0+), or a suitable 
combination of a sufficient number of both ground state and exited 0+ state transitions. In this 
multiple experiment analysis, double electron capture can be come a very useful partner. 
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